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Abstract. - We introduce three measures which quantify the degree to which quantum 
systems possess the robustness exhibited by classical systems when subjected to continuous 
observation. Using these we show that for a fixed environmental interaction the level of ro- 
bustness depends on the measurement strategy, or unravelling, and that no single strategy is 
maximally robust in all ways. 



When describing the observation of individual quantum systems, it is crucial to explicitly 
treat the measurement process using quantum measurement theory. The reverse, however, is 
commonly true for classical systems — it is quite usual for experimentalists both to model 
and gather data about classical systems without any reference to classical measurement theory 
(being the theory of Bayesian statistical inference or nonlinear filtering [1,2]). This dramatic 
contrast is only possible because of certain properties which many classical systems possess. 
Three such properties are apparent if we consider watching the motion of a pendulum: 1. 
When we open our eyes, we obtain the information (the location and velocity of the pendulum) 
almost instantaneously, 2. the information we obtain is not exclusive — many people can 
observe the same system, and will all agree upon the results, and 3. the system is relatively 
unaffected by noise so that if we close our eyes for a moment, we can accurately predict what 
we will see when we open them again {}). We will refer to the degree to which a system 
possesses these properties as its degree of robustness under each property. 

The classical robustness of quantum systems in this sense was the subject of a recent 
study by Dalvit, Dziamarga and Zurek [3]. They concluded that, for a given environmental 
interaction, a single measurement strategy will maximize the robustness for both properties 1 
and 3 above. Here we investigate the degree of robustness of two canonical quantum systems 
for the three properties above (and a fourth), under a wide class of measurement strategies. 
We show that this conclusion is, in fact, not warranted. That is, in general, for a fixed 
environmental interaction, different measurement strategies are required to maximize different 
notions of classical behavior, and that this depends not only on the system in question, but the 

(^)Notc that the existence of chaos in a system does not destroy its robustness: in this case the divergence of 
neighboring trajectories is exponential and thus has a negligible effect on the predictability for short times. 
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dynamical regime. There is thus necessarily a trade-off between different types of robustness. 
While this study is of fundamental interest, it is also of practical importance for the feedback 
control of quantum systems [4-6] , as will be discussed. 

To clarify the above and establish our results, it is useful first to establish some concepts 
and terminology. The continuous observation of a quantum system can always be described 
by treating the interaction of the system in question with an environment consisting of a large 
number of degrees of freedom. As a result of the interaction, the environment continuously ex- 
tracts information about the system, and this information can be read by measuring the state 
of the environment. This has been called quantum filtering [7]. Once one has chosen an inter- 
action, one therefore has the remaining freedom to choose the manner in which to interrogate 
the environment, and different methods lead to qualitatively diff'crcnt kinds of measurements. 
These different measurements, which constitute the measurement strategies discussed above, 
are often referred to as different unravellings of the environmental interaction [8] . 

In what follows, by conditional evolution, or trajectory, we will mean the evolution of 
the system as a result of a particular unravelling (U) and by unconditional evolution we will 
mean the evolution of the system under the same environmental interaction, but without the 
measurement of the environment. The unconditional evolution is simply given by the ensemble 
average {E[-]) over all the possible conditioned evolutions. 

Before introducing quantitative measures of robustness which characterize the speed, non- 
exclusivity and predictability of a system, we briefiy revisit the conclusion of Ref. [3] . Using 
our terminology, they maintained that the unravelling that most rapidly purifies the state 
is identical to the one that will tend to collapse the system into a "pointer state" (^). The 
"pointer states" [10] are those that lose their purity most slowly when the environment is 
unobserved. If true, this would imply that the same unravelling is most robust under our 
properties 1 and 3. Ref. [3] also argued that purity loss and fidelity loss were equivalent for 
the purposes of defining robustness. We show that in general, contrary to both of the above 
conclusions, different unravellings are optimal for these various concepts of robustness. 

Measures of robustness of unravellings - We will consider four quantitative measures of 
robustness. The first, purification time, quantifies the rate at which the measurement provides 
information about the system. To define this we allow the system to evolve to its unconditional 
steady state, and then ask how long it takes, upon switching on the observation, for the 
measurement to increase the purity halfway to its long-time value (which, for an efficient 
measurement is unit purity). Formally, the purification time, Tpur, is the earliest solution of 

= 0, (1) 

where 6 = ^(1-|-Tr[pgg]). Here pss is the unconditional steady state of the system and S^it^, ti) 
is the non-linear, stochastic, completely positive map which takes the system from time t\ to 
t2 during a measurement of efficiency r}. It has the properties 

E [S^{t, 0)] = 0{t, 0) = S^{t, 0) = e^*, (2) 

where C is the Lindblad-type superoperator generating the unconditional evolution. 

The purific;ation time is relevant for the feedback control of quantum systems. Feedback 
control is realized by observing a system and using this information, as it is obtained, to 
control the evolution (usually by modifying the system Hamiltonian) . The effectiveness of 
any feedback algorithm is thus ultimately limited by the extent to which the system state 
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(^)We note that a "pointer state" unravelling similar to this was previously considered in Ref. [9]. 
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is known. The purification time therefore provides an indication of the time required for a 
feedback loop to become effective, when it employs the given measurement strategy. 

Secondly, wc wish to qiiantify the extent to which the existence of multiple observers 
interferes with the ability of each observer to track the system. In order for multiple parties 
to observe the system separately, they must divide up the environment between them. This 
looks, to one of the observers, i, as if she has an inefficient measurement, with efficiency 
coefficient rji, such that rji — 1 [11,12]. We will quantify the robustness of the measurement 
scheme against a division of the environment by asking at what value of efficiency the average 
purity of the long-time observed state is halfway between its value at perfect efficiency, and 
that of no observation at all. We will refer to this as the efficiency threshold, and denote it 
by Vthr', formally this is the value of r] such that 

= 0. (3) 

The smaller rythr is, the more robust a feedback algorithm would be in the face of inefficiency. 

Our third measure of classical robustness, mixing time, quantifies the rate at which the 
system becomes unpredictable — essentially it characterizes the sensitivity of the system to 
noise from its environment. To define the mixing time, we allow the observed trajectory to 
evolve for long enough that the state is pure, and that the unconditioned evolution would 
have reached a steady-state, stop measuring the environment, and consider the decrease in 
the purity of the state as time passes. The mixing time is the time at which the purity falls 
halfway from its initial value (being unity) to the value it would take if the system were 
allowed to evolve back to its steady state (being Tr[pgg]). Formally, the mixing time, Tmix, is 
the earliest solution of 

= e. (4) 

This quantifies the predictability of the evolution while the observer is not looking. 

A related concept, introduced previously by two of us and Vaccaro [13, 14], is the extent 
to which an initially conditioned state remains unchanged during a period of unconditional 
evolution. This is quantified by the survival time, Tsut, the earliest solution of 

^lim E [Tr [SY{t, 0)pss x 0{t + t, t)S^{t, 0)p J ] = 6. (5) 

Note that the stochastic map Si {t, 0) appearing twice in this equation is the same map (i.e. 
has the same noise). We also note that Tsur < Tmix- In general this does not measure the 
noise-sensitivity alone, since any deterministic evolution will contribute to changing the state. 
While the survival time is not motivated by classicality like the other measures of robustness, 
it is, like them, relevant for feedback control. Consider feedback with a time delay r. For a 
feedback algorithm designed for r = 0, the performance should not be greatly degraded as 
long as T < Tsur, as the system will not have strayed significantly during the delay. On the 
other hand, deterministic evolution during the delay could be corrected for in the algorithm 
design. In that case, the performance of the feedback control should not be significantly worse 
than for the case of no delay as long as t < Tmix- 

A particle undergoing quantum Brownian motion - We consider a particle in one dimension 
with position q and momentum p in a viscous environment at temperature T . For our first 
example we will consider all continuous Markovian unravellings [15]. In an optical realization, 
this includes all homodyne and heterodyne detection schemes [15]. The stochastic master 
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equation (SME) describing the conditional evolution of a particle under quantum Brownian 
motion (QBM), for all such unravellings, is 

dp, ^ -{i/2)[p' + {qp + pq) /2,p,]dt + V[c]pcdt 

+ {V^dW^(c-Tr[cpe])Pc + H.c.}. (6) 

Here c = ^jTTq + ip/y/WT and 

V[B]A^ BAB'' - i {B'iBA - AB'' B) . (7) 

The subscript c indicates that the evolution is conditioned on the innovation [1] dW , a stochas- 
tic Wiener increment (to which we will return below). Finally, note that we are using scaled 
units such that the damping rate 7, the particle mass m, and h are all unity. Averaging over 
the noise (or setting the detection efficiency rj to zero) removes the second line in Eq. ((BJ, 
leaving a version of the QBM master equation which is in the Lindblad form [16]. It is the 
standard time-independent QBM master equation [17], with the necessary addition of a term 
which generates position diffusion [18]. This system is also equivalent to the damping of an 
optical cavity mode with an added non-linearity. 

The Wiener increment in the above SME satisfies dWdW* — dt and dW^ — vdt, where 
V = re*'^ £ C with r < 1. It is the value of v which determines the measurement strategy, 
and thus our search for optimally robust unravellings will involve optimization over the disk 
parametrized by r and 0. In an optical realization, r ~ corresponds to heterodyne detection, 
and r = 1 corresponds to homodyne detection of a linear combination of q and p, x = 
^g-i0/2 ^tg*0/2 (0 = corresponds to a measurement of q and (j) = n to one of p). 

Since the unconditional (77 — 0) steady-state of Eq. © possesses a Gaussian Wigner 
function, and the conditional evolution preserves Gaussianity, the four measures of robustness 
with which we are concerned may be written solely in terms of the variances of q and p and 
their covariance Cqp = {qp + pq)/2 — {q){p). In particular, the purity, which is required for 

three of the four measures, is P = l/^4:{VqVp — C^p). Furthermore, these quantities evolve 

deterministically even under conditional stochastic evolution. This is very useful, as it removes 
the need to do stochastic simulations for the ensemble averages in Eqs. - l|Sl)- This enables 
us to numerically optimize robustness over all the measurement strategies (all values of v). 

Optimizing for each form of robustness over the disc of measurement strategies, we find 
that the optimal strategy always lies on the boundary of the disc (r = 1), and thus corresponds 
to some form of homodyne detection. However, the type of homodyne detection required (that 
is, the value of the phase (f)) is highly dependent upon the type of robustness desired. The 
results are summarized in Figure ^ where the optimal phase for each notion of robustness 
is plotted as a function of the temperature of the bath. We see that while the purification 
time and efficiency threshold may be optimized more or less simultaneously, the mixing and 
survival times require quite different unravellings (^). 

The two level atom - The unconditional evolution of our second example, a driven two 
level atom (TLA), is governed by the resonance fluorescence master equation 

^ = -i^[a,,p]+V[^a-]p, (8) 

(^)It is worth noting a subtlety regarding the optimal strategy for the survival time at large temperatures. 
While it is true that setting = 7r results in a homodyne measurement of momentum, it is incorrect to 
conclude from Fig. that the optimal survival time is obtained by a momentum measurement for large T. 
This is because, for values of close to tt, the measured observable is given by a; ~ — —^=p+ (tt — (f>)^^p-q. For 

large T, the </> which maximizes Tsur is such that (tt - (/)) ~ T~^^^. Thus X is, in fact, dominated by position, 
as for the other measures. 
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Fig. 1 - The angle, (p, which gives the optimally robust measurement scheme for each of the four 
forms of robustness, plotted as a function of the dimensionless temperature. 



where fl is the Rabi frequency and 7 is the spontaneous emission rate. Since the electromag- 
netic field is the bath which mediates the measurement, all measurement strategies involve 
observing the atomic radiation. The equation describing a particular unravelling will include 
additional stochastic terms specific to the unravelling. Note that unconditional dynamics of 
the TLA are entirely determined by the dimensionless parameter n/7. 

For the TLA we will consider both continuous unravellings and those which involve discon- 
tinuous jumps. However, in this case we cannot eliminate the need to take ensemble averages 
over large numbers of trajectories (for the results which follow, typically hundreds of thou- 
sands). As we are not therefore able to optimize over all unravellings, we choose a small set 
which are most relevant, either from an experimental point of view, or because they posses 
special properties. We consider direct photon counting, homodyne x detection (0 = 0), ho- 
modyne y detection ((/> = tt), heterodyne detection, and Adaptive Interferometric Detection 
(AID). This last measurement scheme was introduced by Wiseman and Toombes [19], and 
involves interference of the emitted radiation with a local oscillator (LO) as in homodyne 
detection. However, unlike in homodyne detection, the LO is weak (comparable in amplitude 
to the TLA field) so that individual photons are resolvable, resulting in a jump process. Upon 
each jump the amplitude of the LO is flipped via a real-time feedback loop, which makes it 
adaptive and non-Mar kovian. We consider AID because it has been shown in reference [14] 
that it is the optimal unravelling for maximizing the survival time, and is thus a likely candi- 
date for maximizing the other measures of robustness in which we are interested here. 

We now evaluate the three measures of classical robustness for each of the above mea- 
surement schemes, to rank them. The fourth robustness measure, the survival time, has been 
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Table I - Rankings of robustness of unravellings for each measure of robustness with most robust at 
the top down to least robust at the bottom. 
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calculated for this system previously in [14], and we include these results here for comparison. 
The rankings of the unravellings are displayed in table HJ from most to least robust. 

The ranking of the unravellings determined by the survival and mixing times are identical, 
and independent of il/7; AID remains the most robust by these measures for all dynamical 
regimes. The ranking under purification time depends on the regime. For weak driving 
(fi < 7/2), homodyne-x and AID provide the most rapid means of obtaining information 
about the system. However for strong driving AID is one of the least effective at extracting 
information. The full dependence of the purification time for all the schemes is displayed 
in figure [3 The ranking in terms of the efficiency threshold also depends on the dynamical 
regime as can be seen in TableQ] For 7 $7 (not shown in the table) direct detection actually 
gives the most robust efficiency threshold. 

Discussion - We have introduced means of quantifying the degree to which quantum sys- 




D.J. Atkins, Z. Brady, K. Jacobs and H.M. Wiseman: Classical Robustness of Quantum UnravellingsT 



terns behave in a classical fashion under continuous observation. We have also investigated, 
for two canonical systems, how this classicality or robustness depends on the way in which 
the environment is interrogated. 

Reviewing the results for both systems reveals that for the most part, of all the Marko- 
vian measurement schemes, homodyne detection provides the most classically robust means 
of observing the systems. However, for QBM the different concepts of robustness require dif- 
ferent homodyne schemes (i.e. different (/)). For the TLA, an adaptive (hence non-Markovian) 
measurement strategy is most robust in general. Moreover, for the TLA, there are dynamical 
regimes (when looking at purification time and efficiency threshold) where direct and het- 
erodyne detection render the most robust behavior. In summary, it is clear that there is no 
unique unravelling which is the most classically robust, contrary to previous expectations [3]. 

The notions of robustness which wc have considered here also provide an indication of 
the relative merits of different measurement schemes for feedback control constrained by time 
delays and measurement inefficiency. The fact that no single unravelling is maximally robust in 
all ways suggests that the measurement strategy adopted for the purposes of feedback control 
will need to be tailored to the robustness requirements of a given application. Quantifying 
this link between classical robustness and quantum control should provide useful insight into 
the design of quantum feedback algorithms. 

* * * 
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